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Abstract 

Recent series of articles (B.M. Tanygin, 2011-2012) has been aimed to 
study flexomagnetoelectric properties of magnetically ordered crystals by 
means of a magnetic point symmetry description. Besides its fundamental 
interest (complete symmetry classification), applications to a classical micro- 
magnetic variational problem of a thermodynamically equilibrium state cal¬ 
culation had been provided without a proof. It was shown that the magnetic 
point symmetry group determination of the given micromagnetic structure 
could be used to predict qualitatively the electric polarization distribution 
induced by a flexomagnetoelectric coupling. This article is an important 
addendum which provides a formal mathematical explanation why and how 
a symmetry-based description of flexomagnetoelectric phenomena conforms 
to a quantitative calculation of an equilibrium state of a magnetization and 
electrical polarization vectors spatial distribution. An application of the Neu¬ 
mann’s principle and a group-theoretical methodology to real micromagnetic 
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structures is clarified as well. The well-known role of a symmetry-based 

approach in new interactions and phenomena discoveries has been discussed. 
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1. Introduction 


A group-theoretical approach was originally well applied to ferroelectric 


domain structures 


l|. A magnetic point symmetry classification of different 


2-|l3]: the complete list of 


micromagnetic structures had been established 
magnetic point groups is obtained for magnetic domain walls, Bloch lines, 
Bloch points, magnetic vortex and skyrmion structures. An important prac¬ 


tical feature o 


tively 


this symmetry classification is its ability to predict qualita- 


3 a magnetization M and/or electric polarization P vector 
equilibrium spatial distribution: an even/odd/zero/other type of components 
Mi ( Tj ) and P{ (r^) had been identified for each point group. 

Hence, a knowledge of the magnetic point gromyGp of a crystal sample 


simultaneous paramagnetic and paraelectric phase [14] allows to describe (by¬ 


passing calculations) equilibrium spatial distributions of magnetization and 
polarization with 


9] or without 8] their magnetoelectric coupling. Although 
this statement is very natural in the framework of the group-theoretical con¬ 
sideration, it requires a proof in the scope of the micromagnetics 15] and 


the Ginzburg - Landau theory 16]. Also, quantitative parameters cannot be 


obtained from a symmetry theory. The quantitative calculation of an equi- 
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librium state of a magnetization and electrical polarization vectors spatial 
distribution is a solution of micromagnetic problem, which is based on a free 


energy functional minimization 
geneous magnetoelectric coupling 
(FME) coupling 


15, 16]. On the other hand, the inhomo- 

17] (also known as flexomagnetoelectric 


symmetry 


18]) free energy density expression is defined by a crystal 


19], Hence, two approaches are obviously related: a variational 


problem solution and a symmetry-based prediction. This relation was shown 
only for some micromagnetic structures without a general proof. 

The derivation of this proof as well as a methodological analysis is the 
purpose of this work. As a constructive example, only Bloch lines in a planar 
domain wall has been considered. A generalization to other micromagnetic 
structures is straightforward. 


2. Theory 

A calculus of variations of a total free energy functional <f> is given by: 

5 J[F M (r, M, V X M, V y M, V Z M) + 

+Fp (r, P, V X P, V y P, V Z P) + 

+Fme (r, M, P, V X M, V y M, V z M)]d 3 r = 0, (1) 


where XYZ is the reference Cartesian crystallophysical coordinate system 

20]; the ferromagnetic Fm, ferroelectric Fp, and flexomagnetoelectric Tme 
19] components of the local free energy density is integrated over the whole 
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crystal sample volume. The equation (JT]) could be supplemented by addi¬ 
tional conditions, like the equality |M| = Mg for a ferromagnetic material. 


This problem conventionally reduces to a set of Eulcr- 
which can be solved numerically in the general case 


iagrange equations 


21]. The FME inter¬ 


action free energy density (here and hereinafter, Einstein notation is being 
used): 

Fme = lijkiPiMjV kMi ( 2 ) 

is an invariant of a magnetic point group of a spatially unlimited crystal in 


a paramagnetic and paraelectric phase G 


lilli: 


9 ■ % jk iPiMjV k Mi — r j i j k iP i Mj'V k Mi, 


( 3 ) 


Here and hereinafter, a magnetic point symmetry transformation g G Gp* is 
applying to individual vectors and an operator V, which is treated as a polar 
time-even vector. The FME interaction tensor 7 ij k i is a material constant. A 
total free energy minimization in scope of the variational problem ([Tj) leads 


to the following expression 


22 


,23]: 


Pi (r) = Xe (r) 7 ijkiMj (r) V k Mi (r) , (4) 

where y e (r) is a dielectric susceptibility of a crystal paramagnetic phase. 

In case of a spatially restricted crystal sample, the magnetic point group 
Gp of a crystal sample in the simultaneous paramagnetic and paraelectric 
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phase 14] is derived by means of the Curie symmetry principle: 


Gp = Gp n G s , (5) 

where G s is a magnetic point group of a crystal sample shape. In case of a 
thin him or bulk plate, it is given by the limit Curie group supplemented by 
a time-reversal operation 1': 


G g = oo/mmmT 


( 6 ) 


The symmetry transformations g G Uj of a Bloch line magnetic point 
group Uf. 

Ut C G P (7) 


correspond to functional dependencies of order parameters 


, 4, [l3|: 


g • M (r) = M (g • r) 
9 • P (r) = P (g • r) 


( 8 ) 

(9) 


The detailed algorithm and the table data of the magnetic point group deriva¬ 
tion depending on the boundary conditions and the domain wall plane as 
well as Bloch line orientation had been published before 8), |ll|, 14]. Ac¬ 
cording to (J5]) and ((7]- 0, the crystal sample free surfaces symmetry (Gs) 
perturbs the spatial distribution of electric polarization and magnetization 
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vectors. Although a magnetic point symmetry of this perturbation can be 
the same for the same experimental sample shape like (l6lh an ab initio mech¬ 
anism of this impact can be different: the strain-induced magnetocrystalline 
anisotropy, the demagnetization and depolarization phenomenon, the surface 
anisotropy 24], the boundary conditions of inhomogeneous magnetostrictive 


and piezoelectric effects, the non-local magnetoelectric coupling 4], etc. 

As a rule, each dependency M* (rj) and Pi ( rj ) (any combination: i = x,y,z 
and j = x, y, z) is a sum (F) of an even (S) and an odd (A) function (the 


notation 


11 ] is being used) in the triclinic-pedial (Gf? = 1') magnetically 


ordered crystal. In case of an arbitrary crystal, aside of the transformation 
1 ', additional crystallographic magnetic point symmetry transformations g 
apply sequentially restrictions (0 0 to the components of the dependencies 
M (r) and P (r) compare to the lowest symmetry case. Possible transfor¬ 
mations g are n, n', n, and f\!\ where n e {1,2, 3,4, 6}. This leads to the 
dependencies M t (rj) and Pi (rj) symmetry increase: 


(F) (S) 

(F) (A) 


( 10 ) 


Let us consider (here and hereinafter) only classical planar micromagnetic 


structures (scope of 


Ui C mmml' 


ll|). i.e. a planar domain wall with a Bloch line inside: 


13|. Then, the group Ui consists of the transformations g with 


n — 1 or 2 only. Moduli of the vector components Mi and P,; are not changing 
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by such symmetry transformations. For each multi-index ( ijkl ) (each index 


can be x, y, or z in the crystallophysical coordinate system 


20 ]), both sides 


of the equation (0J) changes sign by the transformation g only simultaneously 
due to the requirement (J3]) taking into account (J7|). 

Hence, the equation (J3J) satisfies the components functional dependence 
parity requirements mm- Former one is a predictive capability in scope 
of a magnetic symmetry consideration of the given crystal. Latter one is 
a solution of the FME coupling variational problem ([T]): a magnetization is 
treated as a known function and a polarization is treated as an unknown one, 
which means that a stable micromagnetic structure induces a polarization 
by FME coupling. Important to note, that any requirements satisfaction 
by the equation (|4|) does not mean that a solution of (JTj) cannot have other 
symmetries, which will be considered in the next section. 


3. The Neumann’s principle application 

If some component of M (r) is zero in terms of a total free energy mini¬ 
mization and its components type of a spatial parity cannot be found in the 
tables with the magnetic point groups list {llj] then the Bloch line group Ui 
should be selected in a way taking into account the rule that a zero functional 
dependence is a particular case of an even or odd dependency: (0) —> (H) or 
(0) —> (S) or (0) —> (. F ) for any component M % (■ rj ) and Pi ( rj ). In this case, 
the group Ui is a medium symmetry after a spontaneous symmetry breaking 
Gp Ui through a magnetic ordering; and a magnetic point group V of the 
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variational problem solution [M(r).P(r)] represents a physical properties 
symmetry (c/. Neumann’s principle 25J, [2G]): 


VDUtCGp 


( 11 ) 


Here, the group V describes completely a magnetic point symmetry of equi¬ 
librium magnetization and polarization spatial distributions in a framework 
of a given model (JTj) independently on a crystal sample symmetry Gp. As an 
example of this concept: a magnetic vortices and skyrmion point symmetry 
has been described by a rotational infinity-fold symmetry axis oo z 12] even 


though this symmetry element does not exist in crystallographic groups 


20 ]. 


Same methodology should be followed in the case when a particular set 
of even and odd functions Mi (rA and Pi ( rj ) parity types cannot be found in 
the tables published in the Ref. ll|, 13]. According to Neumann’s principle, 
following symmetry reduction is possible in case of group Ui selection: (A) —» 
(. F ) or (S) —> ( F ) for any component Mi (rj) or P* (r 0 ) in the given Bloch 
line. It means that an even and odd dependency is a particular case of an 
arbitrary dependency ( F ). 

It is important to note, that as it was shown in the Ref. 3], if a variational 
problem solution provides a too highly symmetrical solution then it means 
that a some free energy invariant has been omitted. Hence, described here 
methodology of the magnetic symmetry-based predictions application to the 
real variational problem solutions can be used as a pointer towards the new 













material properties search. 


Let us consider a similar situation in case of a FME coupling. The well 


known expression of the FME coupling 


22 


, 23, [28] is given by a Lifshitz 


invariant-like single-constant coupling term: 


^me — 7 ijkiPi (Mj-VfcMj — M{S7 fcMj), 


( 12 ) 


which can be derived from the equation (J2]) by a removal of total spatial 


derivatives It was shown 


10 ] that this step cannot be justified by 


the medium symmetry invariance requirements ([3]). Hence, a reason of a wide 
usage of the expression (fT2l) is related to an experimental experience; i.e. a 
predominance of a specific microscopical mechanism in case of some published 
experimental observations of the FME phenomena. A microscopic nature of 


an exchange-relativistic 


4] FME interaction has been suggested based on the 


spin supercurrent model with the similar expression for the cases 


constants) of the superexchange and double-exchange interactions 29]: 


different 


P = -r [e 12 x [Si x S 2 ]], 


(13) 


where Y is a material constant; Si and S 2 are spins of transition metal ions 
in the M — O — M three-atom triad. The e 12 is a polar time-even unit vector 
directed along the line segment connected metal ions. 

The described type of the FME coupling (TT2h requires two non-zero mag¬ 
netization components. According to the term (TT2lh magnetic domain walls 
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or spin-density-wave with a collinear magnetic ordering cannot induce an 


electric polarization. Oppositely, according to the 


1 ly . if at least one mag¬ 


netization components M a (r) 7 ^ 0 then nonzero polarization components 
can be induced by the FME coupling by means of the free energy density 
term 19]: 

fl = licxiaPi^ iMl/2 (14) 

in case of Gf? belonging to a cubic, orthorhombic system, or tetragonal sys¬ 
tem (dihedral groups only) 19]. A formation of multiferroic phases at the 
magnetic commensurability transitions of YM 112 O 5 has been explained by the 
provided FME coupling type (1T41) as a phase dislocation related mechanism 
23], which is suppressed by the well-known expression (fl2l) as a total spatial 
derivative term nn Hence, a progress of an experimental and theoretical 


research of a novel FME coupling mechanism leads 
as a symmetry based one (J4j) published originally 


s 

to 

4, 

17 


o the same conclusion 


Generally, according to the Neumann’s principle, a zero value of a specific 
component Pp (r) is possible in the variational problem solution with given 
model assumptions. Example is a saturation magnetization M (r) = const = 
Mg constancy in a case of a ferromagnetic material leading to the suppression 
of the invariant (TT4l) in crystals of a cubic system 19]. Also, it is possible 
in a case of a given measurement tolerance and a macroscopic level of con¬ 
sideration corresponding to a physically infinitesimal volume selection in a 
research program. Here, if a symmetry-based prediction allows Pp (r) 7 ^ 0 
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then a physical mechanism of the FME coupling should be searched: an 
increase of a variational problem solution accuracy; a considering of all non¬ 


zero (allowed by GfP) values of material tensor components 7^; looking for 


,Q 


expressions of a total free energy additionally 


a new integro-differential 
to the ([3]). 

The Bloch line always contains direction of a continuous translational 
symmetry. If g ■ ( e y r ) = — (e^r) and g ■ AP = — Mi or g ■ Pi = —Pi then 


corresponding components should be zero 


points 


11]. This is not the case for Bloch 


13], which are more general micromagnetic structures compare to 


Bloch lines. In other words, a Bloch line is a particular case of a Bloch point. 
The Neumann’s principle application takes place here as well. Magnetic 
point group Ui of the Bloch point never requires M a (r) = 0 or Pp (r) = 0. 
It was shown that that the set of all possible groups Ui (magnetic point 
symmetry classification) of Bloch points is identical to the one of the Bloch 


lines 


13] . There are 48 magnetic point groups of the Bloch points including 


22 (11 time-invariant and 11 time-noninvariant ones) enantiomorphic and 26 
non-enantiomorphic groups. 


4. Conclusion 

Hence, symmetry-based predictions of the FME phenomena correspond 
to a variational problem solution. Neumann’s principle corresponds to re¬ 
moval of negligible interaction mechanisms or a model related coarse-grained 
approximation. It can be treated as a designation for further discovery of new 
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physical interactions in magnetically and electrically ordered media including 
a coupling of these two subsystems. 
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